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The four-state model is proposed for description of phase transition in 
ferroelectric crystals of DMAGaS and DMAAIS type. Thermodynami- 
cal functions of the model are obtained in the mean field approximation. 
The phase transition between paraelectric and ferroelectric phases is in- 
vestigated. It is established that order of the phase transitions depends 
on relations between model parameters. 
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1 Introduction 

A peculiar feature of isomorphous crystals with ferroelectric properties 
(CH 3 ) 2 NH 2 A1(S0 4 )2 • 6H 2 (DMAAIS) and (CH 3 ) 2 NH 2 Ga(S0 4 )2 • 6H 2 (DMAGaS) 
is the possibility of existence in three different phases at change of temperature: at 
room temperature crystal is paraelectric, at lowering of temperature it sequentially 
becomes ferroelectric and further can be in antiferroelectric state |l|, 0, [3|. [| . 

It has been measured T c \ = 136 K and T c2 = 113 K [|J (or 122 K and 114 K respec- 
tively 0) for DMAGaS crystal, but it has been only found T cl = 150 K |[ |, | for 
DMAAIS. The phase transition between the ferroelectric and antiferroelectric phases 
is of the first order |3], |j] . The phase transition paraelectric - ferroelectric is of the first 
order close to the second one in DMAGaS and of the second order in DMAAIS. Crys- 
tallographic analysis shows that in all three phases crystal belongs to monoclinic space 
groups: high-temperature paraelectric phase has P2i/n space group ^ |j, ferroelectric 
and antiferroelectric ones have Pn [fj, |3| and P2i |3j groups respectively. It should be 
mentioned that low-symmetry space groups are subgroups of the high-symmetry group 
obtained by the loss of rotation axis and mirror plane, respectively (point symmetry 
group 2/m changes to m or 2). 

Data of investigation of various physical properties of the crystals were presented 
in a range of publications. First of all it should be mentioned dielectric, pyroelectric 
and dilatometric measurements (|, H ||, |7|, H, H HI HH HH 5 ultrasonic and NMR inves- 



tigations | 13| , |i~4 , [T5H , studies of relaxational dynamics and lattice dynamics |16|, |17]] by 
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Figure 1: Projection of N and C atoms in the dimethylammonium groups of 
DMAGaS crystal at ambient temperature |l[] onto the XZ plane (position in- 
dices are indicated in parenthesis). The Y coordinates of atoms fol- 
lows: Yjv(i,i)=-0.1469 A, Yjv(i,2)=0.1469 A, Y 7V( i i3) =-0.6736 A, Yjv(i, 4 ) =0.6736 A, 



Y JV(2 ,i ) =5.5099 A, 



Y 



7V(2,2) : 



=5.2161 A, Y 



7V(2,3) : 



=4.6894 A, Y 



7V(2,4) : 



=6.0366 A, Y c(1)1) : 



-0.2606 A, Y c( i >2) 



=0.2606 A, Y 



C(2,l) : 



=5.6236 A, Y 



C(2,2) : 



=5.1024 A. 
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Figure 2: Projection of N atoms onto the plane perpendicular to the C-C axis of the 
dimethylammonium group. 
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means of radiospectroscopy and Raman spectroscopy methods. Such issues as dielec- 
tric anomalies in the vicinity of phase transition points (mainly near T cl ), clarification 
of the order of phase transition, dynamics of ionic groups and protons on hydrogen 
bonds in the wide range of temperatures |L6], |T?|| , influence of external hydrostatic 
pressure 0, [19], |(| were considered. Temperatures of phase transitions were justified; 
according to the latest data T c \ = 136 K, T C 2 = 117 K (DMAGaS) and T c \ = 155 K 
(DMAA1S). 

Crystals DMAA1S and DMAGaS belong to the type of ferroelectrics which have 
some structural elements which reorientation leads to polarisation of crystal. In this 
case the element, which can be reoriented, is the dimethylammonium cation (more 
strictly NH2 group). Crystals contain two symmetrically nonequivalent DMA groups 
per unit cell (Fig. [[]). 

DMA can occupy four equilibrium positions which are related in pairs by inversion 
centre forming a slightly deformated cross (Fig. 0). It is very likely that this asymme- 
try (and corresponding difference of energies of interaction between groups in various 
positions) is responsible for so complicated behaviour of the crystal. In the paraelec- 
tric phase the site in one pair is occupied with probability 40% and in another with 
probability 10% at 300 K. 

Possible role of reorientations of DMA groups in phase transitions of the considered 
family of crystals was pointed out in a range of articles p, || [14], ^TJ. This assumption is 



directly confirmed by the data of structural investigations ]T4|. [2l|] . Performed recently 
NMR measurements [[H] have proved that reorientations and orderings of DMA groups 
can be considered as an origin of ferroelectrics (F) and antiferroelectrics (AF) phase 
transitions while other changes in a lattice structure (rotations of SO4 groups, freezing 
of protons in certain positions on hydrogen bonds) are just accompanying phenomena. 

So for theoretical description of phase transitions to ferro- and antiferroelectric 
states in DMAA1S and DMAGaS crystals and development of appropriative micro- 
scopic model one should take into account reorientation of DMA groups on and their 
influence on the state of crystal. Such an approach is a main goal of the present paper. 
Our task is to formulate this type model, to introduce order parameters responsible for 
appearance of low-symmetry phases, to calculate thermodynamical potential proceed- 
ing from the model and to investigate thermodynamically stable states. The transition 
from paraelectric to ferroelectric state and the dependence of the order of this transition 
on the system parameter values will be studied in more detail. 



2 Hamiltonian of the model 

In both high- and low-temperature phases ionic groups belong to two different sublat- 
tices (an elementary cell consist of two translationally nonequivalent groups). We shall 
characterize orientations of the groups by the spatial localization of nitrogen atoms be- 
longing to the groups. Four possible positions of N atom are described by the Hubbard 
projection operators X^,, where n is the number of a lattice site (an elementary cell), 
a = 1,2 is the sublattice index, p = 1, ... ,4 is the position number (X^, = or 1; 
SpA^na = 1). Energies of different positions are pairwise equivalent: S\ = s 2 , £3 = £4 
(see notations in Figs. [I] and ||). 
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The model Hamiltonian can be written down as follows: 

nap nn' eta' 

pp' 

where the first and the second terms describe energy of separate complexes and pair 
interaction J^ a i{nn') depending on their orientational states, respectively. 

Further consideration will be performed in the mean field approximation. By sub- 
stitution XX -> 2(X)X - (X)(X) in (P one can obtain 

H UF = NH C + J2l £ p ~ f* P ] X Z, (2) 

nap 

where the mean field is introduced 

u= E^'W(^) (3) 

n'a'p' 

and 

H < = 2 EEiwra^i) (4) 

n' act 1 
pp' 

Averages (X%£) do not depend on cell number ((X^) = ( X a P )) due to translational 
symmetry which does not change at the transition to low-symmetry phases (e.g. lat- 
tice period is not doubled in antiferroelectric phase). The averages have meaning of 
occupation numbers of corresponding orientational states and can be grouped in linear 
combinations 

^ = ]TtWm (5) 

ap 

which transform according to irreducible representations of point symmetry group 2/m 
of high-temperature phase. Coefficients of transformation (||) are presented in Tab. [1|. 
As one can see all possible irreducible representations (A g , B g , B n , A u ) are realized, 
two linear combinations of a type of (^) belong to each of them. 

Combinations which are transformed according to B u and A u representations, 
are connected with differences of occupations of (a, 1), (a, 2) and (a, 3), (a, 4) positions 
and can play a role of order parameters for phase transitions to ferroelectric (B u ) and 
antiferroelectric (A u ) phases. 

As one can see in Tab. [j] the first of parameters belonging to the B u representation 
(y + ) describes ferroelectric ordering along the ferroelectric axis (the OX' axis in the 
crystallographic plane (ac)) while the second one (w_) corresponds to antiferroelectric 
ordering along OY (the crystallographic axis b). For the A u representation the ordering 
is inverse: the parameter y_ describes antiferroelectric ordering along OX' and the 
parameter u + defines ferroelectric ordering along OY. Thus, the states of ferro- and 
antiferrophase observed in DMAA1S and DMAGaS crystals in reality are mixed. 

Using (|3|) and (0) one can write down 

fap ^ )ap,pj \lp'£,u' i (6) 
a^' 
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where 



,-pp 



otp a'p' 



(7) 



and 



3 Free energy and equations for order parameters 

Proceeding from expressions (|2|), (|6|) and ([|) one can easy derive main thermodynamical 
functions of the model in the mean field approximation and construct a set of equations 
for selfconsistency parameters 

Corresponding partition function is 
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where 



^2 



cosh/3(/Cj /+ + £;,,_) + cosh /?(£;„_ + fc uH 
Q^T cosh (3(k y+ - k y J) + cosh - fc uH 



and 



Qx 


= e~ Pei exp[-(3(k x+ ±k x .)] 


, Qt = 


- e Pez exp[— 


k x + 


= (A lX+ + C lZ+ )/2, 


k x - 


= (A 2 x- 


+ C 2 z-)/2, 


k z + 


= (C lX+ + B lZ+ )/2, 


k z - 


= (C 2 x_ 


+ B 2 z_)/2, 


ky + 


= (A 3 y + + C 3 u_)/2, 


ky- 


= {A 4 y- 


+ C A u + )/2, 




= (C 3 y + + S 3 w-)A 




= {C±y- 


+ fi 4 w+)/2. 



Here such notations are introduced 

A 1 = a + b + g + h, B x = 
A 2 = a + b- g-h, B 2 = 
A 3 = a-b + g-h, B 3 = 
A4 = a — b — g + h, B± = 
a = Jn(0), & = Jn(0), c - 
e = j!K0), / = Ji 3 i 4 (0), 9 



J 



jim, 



Jim, 1 



e + f + l + m, Ci 
e + / — I — m, C 2 = 
e- f -l + m, C 3 = 
e — f + I — m, C4 = 
= Jn 3 (0), d = jl*(p), 

=jim, h=j%(o), 
-jim,™=jim- 



c + d + j + k, 
c + d — j — k, 
c — d — j + k, 
c — d + j — k, 



(9) 
(10) 
(11) 



(12) 



(13) 



(14) 



(15) 
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Table 1: Coefficients of symmetrized occupancies of orientational states which corre- 
spond to irreducible representations of the point symmetry group 2/m. 



Free energy per site F is 

F 



H c -Q\n{2Z l ) -61n(2Z 2 ). 



(16) 



By direct averaging or with the use of free energy ( |TB| ) one can obtain the set of 
equations for the parameters 



x± = [Z^[ l Q+ cosh (3(k y+ + k y J) ± Z 2 X Q~ cosh f3(k y+ — k y J)}/2, 

z± = [Z Q+ cosh f3(k u+ + k u _) ± Z 2 X Q~ cosh/5(/c u+ - k u -)]/2, 

y± = [Z{ l Q+ sinhj3(k y+ + k y J) ± Z 2 X Q~ sinh j3{k y+ - k y -)]/2, 

u± = \Zy X Ql sinh (5(k u+ + k u _) ± Z 2 l Q~ sinh (3{k u+ - k u _)}/2. 



(17) 



4 Phase transition between paraelectric and ferro- 
electric phases 

Let us consider phase transition (PT) from paraelectric phase to the ferroelectric one 
which corresponds to irreducible representation B u with the two-component order pa- 
rameter (OP). In this case averages x_, z_, y- and u + are equal to zero, averages y + 
and U- are the components of OP and Z\ = Z 2 . Then the expression for free energy 
(|l6l) simplifies 



F 



H c -2Qln{2Z l ), 

I (A x Nl + 2CiiV 1 iV 2 + B x Nl + A£ + 2C 3 £ V + B 3 r, 2 ) 



e _/3ei exp 



--PiAxNt + dNz) 



+ e /3e3 exp 



■-Pi&Nr + B,N 2 



cosh 
cosh 



-(3(C^ + B 3V ) 



Here 



iVi = x + ,N 2 = z + ,£ = y+,r] = u- 



(19) 
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It is convenient to separate the term in the expression for the free energy which depends 
on £ and r\ 



F 



Fo + F^v), 

- {A^l + 2C X N X N 2 + B x Nl) - 29 In 2 - 26 ln(a + b), 
(A 3 e + 2C 3 £t7 + B 3 rf) 



(20) 
(21) 



2B In < n x cosh 



-f3(A 3 £+C 3V ) 



+ n 2 cosh 



-/3(C 3 ^+B aV ) 



where the following symbols are introduced 



a = e /3ei exp 



--/3(A X N X + C X N 2 ) 
rix = a / (a + b) , n 2 = b/(a + b), n x 



b = e /3£3 exp 



n 2 . 



--P{C X N X + B 1 N*) 



(22) 



Expressions ni and n 2 describe an effective occupancy of equilibrium positions 1,2 
and 3,4, respectively, in disordered state. One can assume they have a weak dependence 
on temperature in the region of the paraelectric-ferroelectric phase transition. In the 
framework of this approach n x and n 2 could be considered as model parameters. The 
main contribution to the change of free energy is given by the term F x (£, if). Then the 
averages £ and rj are given by a set of equations 



n 



[l-n) sinh[/3(£ + Cr/)] 



n) cosh[/3(£ + Crj)] + n cosh[/?(C£ + Br])} 
nsmh0(C£ + Bri)] 



(23) 



(1 - n) cosh[/?(£ + Crj)] + n cosh[/?(C£ + Br,)} 
where following symbols are used for normalized parameters 



B = B 3 /A 3 , C = C 3 /A 3 , = (3A 3 /2, 



n 



n 2 . 



The set of equations (^) is invariant with respect to the inversion operation £, — > 
— £, — r] and has, as usually, a trivial paraelectric solution £ = 77 = 0. It can be 
established that at = two another solutions (£ = 0, rj = 1 at C < B and £ = 1, 77 = 
at C < 1) exist. 

Landau expansion of the free energy is a convenient tool for investigation of PTs of 
the second order (PT2) and of the first order (PT1) which are close to second order. 
But in this case the use of Landau expansion meets with problems due to a negative 
sign of the coefficient at the sixth order of OP for the some range of parameter values. 
For this reason Landau expansion plays an auxiliary role in our investigations. Besides 
that we perform all calculations of £ and 77 averages, free energy and the temperature 
of PT numerically using the expression ( p0|) and the set (p3|). 

Landau expansion up to fourth order was used for calculation of temperature of 
possible PT2 and boundaries of the regions where PT1 can occur. In this case OP is 
two-component, hence for diagonalization of the quadratic form, which looks like 



F 1 ^,r ] ) = Ue + Vr ] 2 + 2S^r ] , 



(24) 
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U= 2 



l-(l + (C 2 -l)n)| 
1 - (1 + (B - l)n) \ 



V =2 



B-(C 2 + (B 2 - C 2 )n) | 



(25) 



one should make the following transformation 

n\ f U\ f= (cos V -sin^\ 

r 2 J \ V ) \ Bm</3 COS</> J 

cos 2^ = (V - U)/y/{V -U) 2 + 4S 2 , sin2y? = 2S/y/(V - U) 2 + 4S 2 . (27) 



Now the expression (21) is as follows 



*2(ri,r 2 ) = A x r 2 + A 2 r 2 , 

Ai = J [(V + U) - ^{V-U) 2 + AS 2 



:28) 



A, 



2 L 



(v + [/) + V(v - u) 2 + as 2 



Coefficients Ai and A 2 are equal to zero at temperatures 
© 12 = [(B-C 2 )(l-n)n\ I 



1 - (1 - B)n) =F a/(1 — (1 + S)n) z + 4C 2 (1 - n)n 



One can obtain conditions on parameters A3, B% and C3 

A 3 > 0, B 3 > 0, A 3J B 3 > C 3 ; B>0, B> C 2 , 



(29) 



(30) 



taking into account that at © = 00 {(3 = 0) the system is paraelectric (coefficients Ai 
and A 2 are positive). If the conditions (|30D are satisfied, the temperature Gi is always 
higher than the temperature G 2 . Hence it determines the temperature of possible PT2 
and variable 7*1 in the expression ( |28|) becomes the OP. 

The region of parameter values, where PT1 takes place, can be established by the 
criterion of a negative fourth order term proportional to r\ of Landau expansion. This 
term can be obtained from the expression 



-l^{(2k l + 2k 2 C i -{l-n)nC 2 ){i 4 



+ (2/^C 4 + 2k 2 B A 



n)n 



B 2 C 2 ) V 4 



(31) 



+ 2C (Ah + Ak 2 BC 2 - (1 - n)n(B + C 2 )) ^ 
+ 2C (Ak x C 2 + Ak 2 B 3 - (1 - n)nB{B + C 2 )) £rf 

+ (l2(Jfei + k 2 B 2 )C 2 - (1 - n)nC 2 (B 2 + ABC 2 + C 4 )) £V }, 

k x = (1 - n)/12 - (1 - n,) 2 /4, fc 2 = n/12 - n 2 /4, 

using the transformation (|26|). 

It is convenient to make a qualitative analysis of PT2 region in the case C — 0, 
when the free energy ( [24]) has a diagonal form and OP is equal to £ or 77. The equation 
B = 1/n — 1 describes the curve, which separates regions with different OP on the 
B — n diagram. This equation can be obtained by setting equal the temperatures at 



S 



1/3 2/3 

n 



0.1 0.2 0.3 

n 



Figure 3: Diagram of order of the 
paraelectric-ferroelectric phase transi- 
tion at C — 0. Filled regions corre- 
spond to PT1, blank regions - PT2. 



Figure 4: Boundaries of PT1 region 
(inside the curves) at B — 5. Solid line 
corresponds to numerical calculations, 
dashed one - the result of Landau ex- 
pansion analysis. 



which coefficients at £ and 77 in the quadratic form (|24|) are equal to zero. It is found 
that the sign of the coefficient at the fourth order of OP depends on n only. Obtained 
diagram of a PT order is depicted in Fig. [| 

Increase of C leads to the narrowing of the region of n where PT1 is possible 
(Fig. ^|). One can see that the region, determined by the sign change, is smaller 
than obtained by direct numerical calculation. This difference can be explained on 
the example of the case C = 0. The point n' is situated on the OP £ — > T] change 
curve of the diagram in Fig. || On the left from the n' point the OP is £, hence PT2 
on £ should takes place before PT2 in 77. But actually PT1 in 77 takes place before 
them. It is illustrated in Fig. [| where dependences of free energy and parameters £, 77 
on temperature are depicted. It is also possible to find PT2 in £ firstly and PT1 with 
jump of 77 at lower temperature. One should note that at C 7^ the OP r\ is a mixture 
of £ and 77. But at considered here values of parameters the main contribution to the 
jump is given by 77. The diagram in Fig. |3] shows that PT1 with prevailing jump of £ is 
typical for the region of small values of B and large n; with prevailing jump of 77 - for 
large B and small n. Increase of C cause a stronger "mixing" of £ and 77, but there is 
the limit C 2 < B (^) of C value. Thus it seems difficult to find in this approach PT1 
with the prevailing jump of £ (i.e. with the prevailing occupancies of positions 1 and 2) 
at small values of n what is characteristic of DMAGaS crystals at PT from paraelectric 
into ferroelectric phase. At small values of n the model exhibits PT2 (close to PT1) 
with prevailing change of £ parameter (Fig. ||). So, we see that despite the simplicity 
of the model its thermodynamics is very complicated. As an example an exotic case 
with sophisticated dependencies of £ and 77 on temperature is shown in Fig. 0. 
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Figure 5: Dependences of free energy Figure 6: Ibid. B = 0.5, C = 0.7, n = 

and order parameters on temperature 0.05. 
at values of system parameters as B = 
5, C = 0.1, n = 0.12. 
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Figure 7: Ibid. B = 5,C = 2,n = 0.05 
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5 Conclusions 



A microscopic approach based on the allowance for different orientational states of 
DMA ionic groups in crystal lattice is proposed in the present work for description 
of phase transition in DMAA1S and DMAGaS crystals. In the framework of a simple 
four-state model order parameters are introduced and calculations of thermodynam- 
ical characteristics are made. The phase transition from paraelectric to ferroelectric 
state is studied at different relations between parameters describing interactions of 
orientational states of different ionic groups. 

An interesting feature which follows from the symmetry analysis is the strict cor- 
relation between the ordering along the 1-2 position axis (ferroelectric one) and the 
3-4 position axis. Namely, as follows from table [TJ, ferroelectric ordering along the 
1-2 axis is accompanied by the antiferroelectric one along the 3-4 axis and vice versa. 

Ferroelectric state in these crystals is simultaneously improper antiferroelectric one 
along an other axis and antiferroelectric state is accompanied by ferroelectric ordering 
along the mentioned axis. Factors affecting order of the phase transition at T c \ are 
studied in more details on the base of the model. At small values of n this phase 
transition is established to be of the second order close to the first order one (close 
to the tricritical point) under approximation that total occupancies of positions (1,2) 
and (3,4) are assumed to be the same as in disordered phase and temperature indepen- 
dent. This behaviour coincides with the known result for DMAA1S crystal but does 
not match with data observed in DMAGaS crystals. For the last case it should be nec- 
essary abandon the approximation made and explicitly take into account temperature 



dependences of n\ and n 2 values with the use of expressions (|22|) . Such an approach 
should also be adopted for description of the low-temperature transition (at T c2 ) to 
antiferroelectric phase and investigation of conditions of its realization. This can be 
the closest perspective of development of the model. 
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